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Abstract 

We propose a general scheme for dissipatively prepar- 
ing arbitrary pure quantum states on a multipar- 
tite qubit register in a finite number of basic con- 
trol blocks. Our "splitting-subspace" approach relies 
on control resources that are available in a number 
of scalable quantum technologies (complete unitary 
control on the target system, an ancillary resettable 
qubit and controlled-not gates between the target 
and the ancilla), and can be seen as a "quantum- 
controller" implementation of a sequence of classical 
feedback loops. We show how a large degree of flex- 
ibility exists in engineering the required conditional 
operations, and make explicit contact with a stabi- 
lization protocol used for dissipative quantum state 
preparation and entanglement generation in recent 
experiments with trapped ions. 
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1 Introduction and Notation 

State preparation problems are of vital importance 
across quantum information processing applications, 
ranging from initialization of quantum computation 
and simulation algorithms in network-based as well 
as cluster-state architectures to use within enhanced 
quantum metrology protocols |12) . If the prepara- 
tion is to be achieved irrespective of the initial state 
of the system, from a control standpoint this trans- 
lates naturally into a stabilization problem. Phys- 
ically, one is compelled to make the quantum evo- 
lution irreversible, by introducing open-system fea- 
tures in either an open-loop or closed-loop fashion 
[2]. Closed- loop control based on quantum measure- 
ments and feedback techniques provides, in particu- 
lar, a very natural and powerful toolbox |20j . with 
feedback from a single ancilla qubit together with 
fast, complete unitary control allowing in principle 
to engineer arbitrary open-system dynamics |9j. 

Quantum state stabilization problems have been 
studied in depth for continuous-time dynamical mod- 
els, for either Markovian "output feedback" (X5J E3 
ITS] or strategies based on state reconstruction by 
quantum filtering [T7J [TU] . It is worth remarking that 
in such continuous-time scenarios the target state 
can typically be reached only asymptotically, in the 
(formal) limit of infinite evolution time. Recently, a 
linear-algebraic framework for analysis and synthesis 
has also been developed for the discrete-time case. 
In particular, for a given indirect quantum measure- 
ment and complete unitary control over the target 



system, it has been shown that pure states are gener- 
ically stabilizable [3]. This is always true for a projec- 
tive measurement, and if the latter is associated to 
a non-degenerate observable, it further follows that 
the desired state can be prepared in a single step of 
measurement-plus-control [5 . 

From an implementation point of view, a prob- 
lem associated with measurement-based feedback 
schemes is that the required control resources need 
not be readily available for many state-of-the-art ex- 
perimental quantum devices. In this work, we in- 
vestigate how to prepare a desired pure state in fi- 
nite time by means of a reasonable set of resources 
for a multipartite qubit system. In particular, while 
still assuming complete unitary control, we shall ef- 
fectively "encode" the whole feedback loop in a finite 
sequence of coherent feedback actions pj \W\ , where 
the controller is itself a quantum system (a qubit) 
and no measurement is involved. For certain exper- 
imental settings, most notably trapped ions, these 
control resources are not only achievable in principle, 
but have been already experimentally demonstrated 
up to 5 qubits [5] . In particular, we will illustrate how 
the "stabilizer pumping" strategy proposed in [3l [11] 
fits into our general "splitting-subspace" framework, 
and how control actions achieving stabilization of ar- 
bitrary pure target states can be explicitly synthe- 
sized. 

We begin by recalling some basic concepts and no- 
tations. Let S be a finite- dimensional quantum sys- 
tem of interest, with associated Hilbert space T~L ~ 
C d . Vectors and linear functionals on % are denoted 
using Dirac's notaion with \ip) and (</>|, respectively 
[15] . Observable quantities on S are associated to 
self-adjoint operators on H, here represented by Her- 
mitian matrices X = X* e $)(H). In particular, the 
state of S is in general described by a density oper- 
ator p G D(H) = {p G Sj(H)\p > 0, tr(» = 1}, 
with pure states corresponding to the extreme point 
of the (convex) set Unitary matrices are de- 

noted by U G U(Ji). The (real) spectrum of an ob- 
servable X represents the set of the possible outcomes 
in the simplest case of a so-called projective (or von 
Neumann's) quantum measurement on S. Suppose 
that X admits a spectral decomposition of the form 
X = J2i XiHi, in terms of a complete set of orthogo- 



nal projectors {IT} on %. According to the basic pos- 
tulates for von Neumann's measurements, the prob- 
ability of obtaining Xi given the pre-measurement 
state p is Pi = tr(IIjp) = tr(IT/TT). Conditionally 
upon the measurement outcome Xi being recorded, 
the (normalized) post-measurement state of S then 
becomes p\i= ^-ITpIT. 

If S consists of multiple (distinguishable) subsys- 
tems Sj, j — 1, . . . ,N, each associated to a Hilbert 
space Hj, the corresponding mathematical descrip- 
tion is carried out in the tensor product space, H — 
Hi ® %2 • • • ® H.N [H], and observables and den- 
sity operators remain associated with Hermitian and 
positive-semidefinite, trace-one operators on "H, re- 
spectively. Given a state p G D(H) of S, the re- 
duced state of one (or a subset) of subsystems may 
be uniquely determined by taking the partial trace 
over the remaining subsystcm(s), for instance in the 
simplest bipartite setting, % = Hi ® we shall 
denote by pi — tr-^ 2 (p) the reduced density operator 
describing the first subsystem alone. 

If information about a quantum system is gathered 
indirectly, through measurements of a correlated aux- 
iliary system, the formalism of von Neumann's pro- 
jective measurements is overly restrictive in general 
and a description of the effect of the measurement on 
the system of interest is provided by so-called general- 
ized measurements. To a set of possible measurement 
outcomes labeled by k, we associate a set of measure- 
ment operators {M^} on the system of interest, in 
such a way that Y^k^-l^-k = I, with / being the 
identity operator. The probability of obtaining the 
fc-th outcome is then computed as pt = tr(M^Mf.p), 
p denoting now the reduced state of the system of 
interest which, after the outcome is recorded, is up- 
dated to p\ k = j- k M kP Ml 

By taking the average over the possible outcomes, 
we obtain a Trace-Preserving Completely Positive 
(TPCP) linear transformation of the state in the form 
of a so-called Kraus map [B], that is, 

p^E(p)=Y J M kP Ml 

k 

The standard case of projective measurements on S 
is formally recovered by choosing = IT-. 
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2 Discrete-Time Feedback Sta- 
bilization 

Suppose that a generalized measurement operation 
can be performed on the target system at times 
t = 1,2,..., resulting in an open-system, discrete- 
time dynamics described by a Kraus map with oper- 
ators {ATfc}. Suppose, in addition, that we are able to 
enact arbitrary unitary control actions on the state of 
S, that is , p i y UpU\ for arbitrary U £li(H), with 
unitary manipulations that are fast with respect to 
the measurement time scale, so that the measurement 
and the control effectively act in distinct "time slots" . 

We can then in principle implement a Markovian 
feedback law, consisting in a map from the set of 
measurement outcomes to the set of unitary matrices, 

U{k) :k^U k EU{U). 

If the above measurement-control loop is iterated and 
we average over the measurement results at each step, 
the net result is a different TPCP map, which de- 
scribes the evolution of the state immediately after 
each application of the controls: 

p(t + l)=J2UkM k p(t)Mlul, 

k 

Controllability and stabilizability for the resulting 
class of discrete-time, closed-loop dynamics have 
been studied in detail in [51 [51 HI [T] . In particular, 
from the results of 4, Ij, it is immediate to see that 
if the following control resources are available: 

(fl) Arbitrary unitary control actions {E/fc}; 

(f2) A non- degenerate projective measurement, as- 
sociated to a resolution of the identity on S, 

{n fc = \<f> k )(MYLi, 

then the system can be prepared in any desired pure 
state in one step. Let pd = \ip)(ip\ denote the tar- 
get pure state: the desired preparation is then sim- 
ply achieved by choosing control operations Uk such 
that Uk\4>k) = l^)- I R control-theoretic terms, this ef- 
fectively implements a quantum dead-beat controller, 
reaching the desired state not just asymptotically but 



in one step. The above is indeed an abstract descrip- 
tion of the most straightforward procedure to prepare 
a given state: first, measure the system projecting it 
onto some known state, and then enact some open- 
loop, controlled transition to steer it to the desired 
state. Despite its conceptual simplicity, this strategy 
may become challenging in practical control scenar- 
ios where the required measurement procedures are 
unavailable (for instance, measurement may be de- 
structive, too slow and/or inaccurate, or not having 
the needed resolution). In what follows, by focusing 
on multi-qubit systems, we show how to achieve the 
same pure-state preparation by replacing the above 
full resolution measurement with the ability of using 
an extra qubit as a fully coherent resettable quantum 
controller [3 [S]. 

3 The Splitting-Subspace Ap- 
proach 

Consider an TV-qubit register, with associated Hilbert 
space Hq = QjHj ~ C 2 , and with {\<f>j)} 2 j=1 de- 
noting the standard (computational) basis of Hq, 
|0i> = I0...00), |0 2 ) = |O...O1),|0 3 ) = I0...10), 
and so on. 

Assume that the following control resources are 
available: 

(si) Arbitrary unitary control actions {U} on the N 
target qubits; 

(s2) An auxiliary control qubit, with Hilbert space 
"H c , that can be reset to a known pure state, 
say |1); 

(s3) Controlled-not unitaries C\ n , C ou t € U (H c £g> Hj) 
between the control and one of the target qubit. 
Without loss of generality, we can pick the first 
qubit (j — 1) and write C out = 1% ® |0)(0| +a x (g) 
|1)(1|, C in = |1)(1|®/ 3 + |0)<0|(8ctx, respectively. 

We next show that a simple approach can be de- 
veloped in order to design quantum circuits able to 
prepare the target pure state pd = \ip){ip\ in a finite 
number of iterations. The first step is to provide a 
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characterization of the target state in terms of a fam- 
ily of splitting subspaces. 

Lemma 3.1 (Splitting Subspaces) Any N-qubit 
pure state \tp) € T-Lq can be described as the unique 
unit-norm vector in the intersection of N subspaces 
Sk of dimension 2 N_1 , that is, 



N 



span{\ip)} = P| S k 



(1) 



k=l 



Proof. It suffices to provide an explicit way to con- 
struct the Sk, see also Figure [T] for illustration. Start 
by relabeling \tpi) := and complete it with 2 N — 1 
vectors so that {\ipk)}k=i is an orthonormal basi^j 
Next, define Si = span{|?/;fe), k = 1,. . . , 2 N ~ 1 }, S 2 — 
span^fe), k = 1, . . . , 2 N -\ + 1, . . . , 2 N ~ l + 

2 N ~ 2 }, and so on for the remaining subspaces. For- 
mally, by defining the matrices 



h 



1 o 
o 1 

s(fc-i) 



1 




in the |^fc)-basis, we can identify the /c-th splitting 
subspace as 

S k = range(IlW). 

By construction, the only vector in the intersection 
of all such subspaces is l^i)- 1=1 

For the important class of stabilizer states |12j . a 
natural choice of splitting subspaces is provided by 
the one-dimensional joint eigenspaces of the associ- 
ated stabilizer operators. However, as showed in the 
above Lemma, there is no need to restrict to this sub- 
set of pure states. 

We also remark that in control-theoretic terms, re- 
quirement (si) corresponds to complete controllabil- 
ity of the iV-qubit in finite time at the level of its 
Hamiltonian description [SJ. In the context of quan- 
tum computation, this is equivalent to requiring ac- 
cess to a (continuous) set of gates achieving exact 

1 Clearly, there is a lot of freedom in doing so, effectively 
equivalent to choosing an element of SU (2 N — 1). This freedom 
can be exploited in reducing the complexity of the unitary 
operations in what follows. 
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Figure 1: Pictorial representation of the constructive 
argument used in Lemma 3.1 to obtain a splitting- 
subspace description for the target state \ipi) = \ip). 



universality |12j . In practice, if unitary operations 
are built out of a discrete set of universal gates, an 
arbitrary target gate U may be implemented by a 
finite-length quantum circuit only within a non-zero 
accuracy. In this case, the finite accuracy of the uni- 
tary control action will carry over to the stabilized 
state, with a worst-case error growing linearly in the 
number of steps. 



3.1 Stabilization protocols 

Building on the above result, the next step is to show 
that the state of the target qubits can be prepared 
in any splitting subspace S in one iteration given 
the control resources specified in (sl)-(s3). If we 
could perform a projective measurement of ITs,IL;, 
the orthogonal projection onto S and its orthogonal 
complement respectively, it would suffice to use a 
feedback law that does nothing if is measured, 
whereas for the other outcome it applies a U± obey- 
ing 



(2) 



Notice, again, that such a control action would be 
highly not unique, since there is at least enough free- 
dom as associated to the unitary mapping from S 1 - 
to S, which is an element of SU(2 N ^ 1 ). The average 
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total evolution would then be: above steps (l)-(4) can be realized in one time unit, 

this yields: 

p(t + 1) = £ s [p(t)} = H s p{t)H s + U ± -n£p(t)Tl%ul, 

iUii , r u ® p(t+i) = tT Hc {u tot [\i){i\® P (t)}ui t } 

and it is easy to see that the support ol p[t + 1) is L ' 

contained in S, as desired: = Usp(t)ILs + U±Ug p(t)Hgll]_ = £s(pt), 



tr[n s p(i + 1) 



= tr[IL s p(t)} 
= tt[U s p(t)} 



tv[u[n s u x iijp(t)iij: 

tr[IL>p(i)] = 1, 



where we have taken advantage of Eq. pL 

Fortunately, the TPCP map £$ in Eq. (|3j) can also 
be implemented without requiring projective mea- 
surements, by using the extra qubit as a quantum 
controller. Let {IV'fc)} denote an orthonormal basis 
for Hq, such that the first 2 7V_1 elements are a basis 
for S. Then proceed as follows: 



(1) 



Perform a unitary on the target qubits, where 
U^lffik) = IV'fc)- I n this way, the state of (in 
particular) the first qubit contains information 
on whether or not the state of the system is in 
S. 



(2) Perform C ou t on the control and the first qubits, 
thereby mapping the information of the first 
qubit on the state of the control qubit; 

(3) Perform C[ n on the control and the first qubits, 
thereby changing the state of the first qubit de- 
pending on the projection of the control qubit 
on|0); 

(4) Perform the change of basis to return to the 
original basis and (if needed for successive steps) 
reset the control qubit back to 

The net effect of operations (l)-(4) is the following 
unitary transformation on % c ® %q: 

Utot — Ujf,Ci n C ou tU^ 

= U^CtoUlUtCautU] 



1 4> in <-> ^ V ip ^ out L/ ^ 

(|l)(l|«g)7 2 iv + |0)(0|<8>t/x)(J2®Hs 



where again Uj_ satisfies Eq. Thus, by tracing 

over the control degrees of freedom and assuming that 



as claimed. The whole sequence of the 4 operations 
(l)-(4) above will be referred in the following as a 
control step, associated to a coherent implementation 
of a single feedback loop in discrete-time. If we have 
N splitting subspaces constructed as in the proof of 
Lemma 3.1 the basis {IV'fc)} to be used at the be- 



ginning of the above protocol can be chosen at each 
control step to be a reordering of the first one. With 
this choice, the desired state-preparation result can 
then be easily established. 

Theorem 3.1 Any N-qubit pure state pd = 
can be deterministically prepared using control re- 
sources (sl)-(sS) in N steps, irrespective of the initial 
state p. 

Proof. Let {Si } be a splitting-subspace description 
for \ip) as in Eq. 0, with £$ t the associated TPCP 
map given in ([3| . Now implement the map £s N ° • • • 
£s 1 (p). By construction, each £$ e maps ^Hil— l ^ 

Sf- onto (n^iSfc) n5<. (rifeLi^fc) = span{|?A)}. 
Thus, at the iV-th control step, any initial state is 



driven into ( HfcLi ^fe ) — s P an {|' ( /')} 



□ 



We stress that, while the controlled- not operations 
entangling the control and the first qubits, as well 
as the resetting of the control qubit, are the same 
at each step (thus in real- world implementations will 
each take a fixed time to be enacted), the U^p will in 
general differ at different steps, since will include 
a reshuffling of the basis elements in order to obtain 
the "nested" subspace structure used for the proof. 
Since, in general, one can see that up to (N — i)/2 
qubit swaps are needed at step £, this can lengthen 
the total implementation time, making the initial 
choice of the (non-unique) basis {IV'/c)} important for 
efficient implementation. 
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3.2 Effective entangling operations 

In practice, engineering exactly the controlled-not 
gates Ci n and C ou t we have assumed may be hard 
to achieve. Luckily, as it turns out, the way of realiz- 
ing the "feedback" map ^ with quantum controllers 
is highly not unique. Even if the controlled-not gates 
are achievable in principle, different choices may be 
more convenient to implement, motivating us to fur- 
ther characterize the degree of flexibility intrinsic to 
our protocol. More precisely, it is easy to verify by 
direct calculation that any pair of entangling gates 
of the following form can achieve the desired task up 
to some additional unitary operation on the target 
qubits: 

C in = (I a ® Vi) (|1)(1| ® h + |0)(0| ® U±) W\h ® %), 

c out = (/ a ® u x )w (d c ® n + o c ® n^) (i 2 ® f/ 2 ), (4) 

where D c , (O c ) are unitary and diagonal (off- 
diagonal) in the standard basis of W c , Uj and Vj are 
unitary operators on Hq and W is an arbitrary uni- 
tary operator on the whole system, respectively. IIo 
(111) indicates the orthogonal projections on the sys- 
tem subspaces corresponding to the first qubit be- 
ing (1). By invoking assumption (si), Uj and Vj 
can be compensated by suitable unitary control: if, 
in the steps (2) and (3) of the protocol described 
above, we now apply (I 2 ® U\)C oxli (I-2 ® U\) and 
(I2 ® Vi)C[ n (l2 €5 V 2 J), respectively, it follows that 



are obeyed. For instance, if Hs,k is the projector on 
the fc-th splitting subspace Sk, and it holds that: 



= (|i)(i| ® i 2N + |o)(o| ® u x ) (d c ® n 5 + o, 



The net reduced dynamics on the system is not af- 
fected by the phases introduced on the ancilla by 
D c , O c , hence tracing out the ancilla we recover ^ 
again. 

In general, failing to compensate the action of 
Uj , Vj as described could potentially slow down (or 
even prevent) the desired state preparation since 

it need no longer be true that fHfc=i ^ i s 
mapped onto (01=1 ^fe ) ^ Sg, a key step in the con- 



Vi,Hs,k = 0, 

ViiLs >fc = n 5 , fe , 



(5) 
(6) 



then there is no need to compensate for V\. In fact, 
([5]) ensures Sk is still an invariant stable subspace for 
the current control step, and (|6| that V acts as the 
identity on it, without compromising the result of the 
previous steps. A concrete example will be provided 
in the following section. Similar conditions can be 
worked out for Ui, Ui,V<z. 



4 Case Study: Trapped Ions 

4.1 Experimental Bell-state pumping 

Recently, a toolbox for engineering the open-system 
dynamics of up to five qubits has been experimentally 
demonstrated using 40 Ca ions trapped in a linear ge- 
ometry [J;. The control resources specified in (sl)- 
(s3) are available through a combination of single- 
and multi-qubit gates ensuring universal scalable ion- 
trap quantum computation, plus access to a control- 
lable dissipative mechanism on an extra ancilla qubit, 
realized via a combination of optical pumping and 
spontaneous emission. Here, we revisit the imple- 
mented two-qubit "Bell-state pumping" protocol in 
the light of our splitting-subspace approach. Exper- 
imentally, a similar approach has been successfully 
employed to create GHZ states on up to 4 qubits. 
Let us denote the four Bell states as 



-^(|00>±|11», |*±) = -1(|01) ±|10». 



vergence proof of Theorem 3.1 Such compensation 
may, however, be unnecessary if additional conditions 



(7) 

The system, initially in an unknown state specified by 
a density operator p, is deterministically prepared in 
the Bell state l^ - ) by realizing in two steps a quan- 
tum operation p H> |^ r_ )(^~|. It is well known that 
Bell states are simple examples of stabilizer state^] 
each of the four Bell states in |7]) may be uniquely 

2 Following standard notation [12] , we shall denote by 
Xk , Vfc, Zk the multi-qubit operators that act as Pauli matrices 
o- x ,<j y , a z on the fc-th qubit, and as identity on the rest. 
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characterized as a joint eigenstate (with eigenvalues 
±1) of a mutually commuting set of two stabilizer 
generators, for instance Y1Y2 and XiX 2 - The con- 
sidered strategy engineers two maps under which the 
system qubit state is transferred from the +1 into the 
— 1 eigenspace of YiY 2 and XiX 2 - 

In order to implement the first map, three unitary 
operations and a dissipative one have been used. All 
act on the two system qubits, along with the ancil- 
lary control qubit, denoted as before with the sub- 
script c. A key role is played by so-called M0lmer- 
S0rensen (MS) entangling gates [H], of the form 
U X 2(6) = exp(-ifS£) and U Y *(9) = exp(-^5 y 2 ), 
where S x and S y denote collective spin operators, 
Sx = Hi X i an d Sy = J2i Y i- Tne mapping steps 
of the experimental pumping protocols are then as 
follows: 

(i) Information about whether the system is in the 
+1 or —1 eigenspace of X X X 2 is mapped by the 
MS gate U X 2 (tt/2) on H c ® Hq onto the logical 
states |0) and |1) of the ancilla (initially in |1)). 



(ii) 



A controlled gate performs a conversion from the 
+1 eigenvalue of the stabilizer XiX 2 to —1 by 
acting on the ancilla and the first system qubitsj^j 



Gin — 



Zi + llXlje® J. 



(iii) The MS gate U X 2(%/2) is re-applied, in order to 
move the state back to the initial basis represen- 
tation. 

(iv) The ancilla qubit is dissipatively reset to state 
|1>. 

Next, the above cycle is repeated, this time us- 
ing [/y2 (71-/2) gates in steps (i) and (iii), with the 
controlled-gate C- m remaining unchanged. 



3 In the experiment, the following sequence is employed: 
C in = C(p) = C/z 1 (a)t/y(7r/2)C/^ 1) (-a)C/y(-7r/2), with 
U^i-a) = exp(i(a/2)X c Xi), U Zl (a) = e iaZ i , and p = 
sin 2 (a). 



4.2 Bell-state pumping and splitting 
subspaces 

In order to show that the above stabilization scheme 
can be seen as an instance of our splitting-subspace 
approach, it is necessary to take a closer look at the 
structure of the MS gate. The MS entangling gate 
is based on pairwise interaction terms and for the 
present discussion it suffices to consider the explicit 
form already given earlieij^J that is, 

(8) 

where the sum defining the collective spin operators 
S x is understood to be performed over all the ions 
involved in the gate. In the experimental Bell-state 
pumping protocol we have just examined, MS gates 
with a phase angle 9 — it/2 are required. In this case, 
the Ux 2 {^ /2) operator in Q can be decomposed in a 
more explicit form with respect to the relevant H. c ® 
Hq tensor-product decomposition: 

u X 2 (|) - u' x (i ® n_i + x <g> n +1 ), 

where U' x is a unitary 8x8 matrix of the form: 



U' x = h 



a -«r/8 



-1 -1 

-1 -1 

1 -1 

-1 1 



u>>. 



The action of the unitary operator U x on the Bell 
states basis Q is described as: 



u x \*-) 



1 



V2 
1 

1 



= 



V2 
1 

7T 



i7r / 8 |$-) 
i7r / 8 |$+) 



4 In its most general form, the MS gate can be 
parametrized by two angles 9 and (f>, (7ms(^> < /') = 
cxp f— ij(cos<f>S x + sin^Sy) 2 ^. 
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where now the unitary operator U Y reads: 



Figure 2: Simplified version (with a single MS gate 
being employed in each control cycle) of the Bell- 
state pumping protocol experimentally implemented 
in trapped ions by Barreiro et al. [3]. 

Hence, U x is not harmful for our purpose since 
it swaps the Bell states in the +1 eigenspace of 
X1X2 and does not change (except for an irrel- 
evant constant), the other Bell states in the —1 
eigenspace. Explicitly, in the Bell basis £>b c ii ; — 
{|$ + ), I*-), I*")}, U' x can be written as: 



tt" 

u X.Bcll — ~^= e 



1 



I7T/8 


' -X 


o 2 ] 




02 


h 



where the symbol O2 denotes a 2 x 2 matrix of zeros. 
We have thus obtained a decomposition of the MS 
gate £/x 2 (7r/2) in a form that includes two terms: 

1. The conditional gate C* ut = 1(g) II_i +X<S)H+i, 
that coherently transfers to the ancilla qubit the 
information on which of the two subspaces the 
system's state is in; 

2. An additional unitary U' x that is not harmful for 
stabilization purposes, since it commutes with 
the projector onto the eigenspaces of X1X2. 

This means that the Ux^{^/2) gate can be decom- 
posed precisely in the form for C ou t given in Eq. Q , 
and is thus a viable operation for implementing a 
splitting-subspace approach to stabilize the desired 
subspace. 

Similarly, the MS gate U Y 2(n/2) can be decom- 
posed in the form: 



U Y 2 



= u Y (z®u , ^ 1 + x®n' +1 ), 



-17T/8 



1 1 

-1 1 

1 1 

-1 1 



Z® U Y . 



Thus, we have likewise found a decomposition of the 
MS gate U Y 2{tt/2) in the form of Eq. Q, involving 
a conditional unitary C^ ut = Z ® n'_ 1 + X ® n' + i 
and an additional unitary U Y with the same roles as 
for the X1X2 case. This means that also U Y 2(ir/2) is 
an admissible entangling operation for the splitting- 
subspace method. 

Note that in the implementation of the experimen- 
tal protocol described in Sec. 4.1 the second appli- 



cation of the MS gates in step (iii) has essentially the 
effect of canceling the action of U' X ,U Y on the tar- 
get qubits. However, it is easy to see that for the 
first map, there is no actual need for this second MS 
gate. In fact, the net effect of the "residual" U' x is 
to swap the states that were in the +1 eigenspace 
of X1X2, but they end up being pumped in the cor- 
rect subspace anyway. Similarly, one can check by 
direct calculation that the action of U' Y on the sub- 
space that is prepared by the first map, namely the 
one generated by {l^ - ), I* - )}, is: 

Uy\$T) = -4=e _W8 |* + ), 



V2 
1 



The desired state l^ - ) is not perturbed: accordingly, 
if we now apply the controlled-cr z to the other state 
we obtain: 

^1*+) = !*-), 

and the desired state preparation is achieved with- 
out the need for step (iii) or, equivalently, without 
undoing U X ,U Y - This simplified version of the sta- 
bilization protocol, which in fact corresponds to the 
experimentally implemented one as noted in the Sup- 
plementary Material [3], is depicted in Fig. [2j 
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5 Conclusions 

We presented a general framework to design sta- 
bilizing controls that prepare arbitrary pure states 
on qubit registers in finite time: once a represen- 
tation of the target state in terms of splitting sub- 
spaces is chosen, the control objective can be achieved 
either by measurements and discrete-time classical 
feedback @J E], or by resorting to coherent feedback 
via quantum controllers ;7 5 . In the latter case, the 
measurement and feedback steps are effectively re- 
placed by suitable conditional operations (C out , C ln 
in our setting), which enable the required quantum- 
information flows out of the system to the controller, 
and viceversa. This fully coherent implementation 
can be practically advantageous in a number of situ- 
ations where quantum measurements are exceedingly 
slow and/or inaccurate, or destructive for the system 
itself. 

The proposed splitting-subspace approach leaves, 
in its current form, significant freedom in construct- 
ing a suitable basis for the system state space. While 
we have shown how the latter translates in added 
flexibility for implementing the required conditional 
operations and can be exploited in principle to min- 
imize the complexity of the stabilization protocol, 
systematic techniques for optimizing the generation 
of a splitting-subspace decomposition in specific con- 
trol settings remain an interesting problem for future 
studies. We have illustrated our general approach by 
revisiting a recently proposed protocol for dissipative 
Bell-state preparation in ion traps [3] in the light of 
an explicit splitting-subspace analysis, demonstrat- 
ing its direct applicability to current experimental 
situations. We expect that this will pave the way 
to further applications of our ideas to synthesizing 
finite-time dissipative state-preparation protocols in 
other qubit devices and control scenarios of experi- 
mental relevance. 
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